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The effect of a thermal gradient parallel or antiparallel to the the main flow direction 
of a pressure-driven symmetric electrolyte in a slit-microchannel is investigated, specifi- 
cally with respect to the elect rokinetic streaming potential. Based on the non-isothermal 
Nernst-Planck equations as well as the Poisson equation, and under the assumption that 
the thermal mobility is the same for each ion species, a correlation of the electric dou- 
ble layer (EDL) potential is derived. The result indicates that the local EDL thickness 
increases (approximately) exponentially with temperature. Consequently, a temperature 
gradient applied along a channel with a constant wall ^-potential leads to a corresponding 
gradient in the EDL thickness. This modifies especially the ion transport by advective 
streaming compared to the isothermal case. Within the Debye-Hiickel approximation, the 
Navier-Stokes equation with the corresponding electric body force term is solved. Ana- 
lytical expressions for the flow profile and the streaming potential under non-isothermal 
conditions are derived. For higher ^-potentials, the results are compared with full numer- 
ical simulations. It is shown that for large pressure differences Apo, the non-isothermal 
streaming potential can be adequately described by the well known isothermal expres- 
sion if the local modification of the Debye length due to the thermal effect is taken into 
account. For a small intrinsic Soret coefficient S* , the total streaming potential becomes 
linearly dependent on the temperature difference AT. For small channel heights at small 
driving pressure differences, the streaming potential is over-predicted (under-predicted) 
by the isothermal expression for positive (negative) values of S'*Ar/Apo- For vanishing 
pressure differences, a nonzero flow velocity and a steady-state thermoelectric potential 
are found. Under such conditions (supplemented by the assumption of local charge neu- 
trality), thermoelectric potentials in bulk electrolytes are known to vanish. It is shown 
that the axial flow and thermoelectric potential found herein results from a confinement 
effect not previously discussed in the literature. 
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1. Introduction 

Over the last couple of decades, electrokinetic fiow phenomena have received significant 
attention by the scientific community. Within the general framework of electrohydrody- 
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namics (Castellanos 19981, the motion of fluids carrying submerged electric charges in 
an electric field as well as the transport of these charges relative to the suspending fluid 
is considered. Electrokinetics is of crucial importance in the stabilization and motion of 
particles in colloidal suspensions (Russell fc Saville|1989 ) and is relevant in electrospray- 
based fabrication methods ( Salataj|2005J or DNA-manipulation/separation techniques 
( ViovypOOO I , to name a few. Ion transport in dilute electrolytes is commonly captured 
with the Nernst-Planck equation. At moderate to high ion concentrations though, the 



finite size and the discrete nature of the ions have to be taken into account (Nadler 
et a/.|[2004 |. To capture the momentum transfer between ions and the solvent, the con- 



ventional stress tensor in the Navier-Stokes equation is supplemented by the Maxwell 
stresses. In case of electrokinetic flow with a free surface, these additional stresses are 
responsible for a number of electric field-induced instabilities of the gas-liquid interface 
as, for instance, observed in the breakup of liquid electrolyte threads subjected to an 
orthogonal electric field (Conroy et al. 2010). Being a manifestation of the Onsager re- 
ciprocal principle (Onsager 1931), the interaction between the suspended ions and the 
liquid carrier goes along with two distinct types of electrokinetic coupling with single- 
phase fluids: one where an electric field drives a fluid motion such as in electro-osmotic 
flow (EOF), and another where ions advected along with the fluid generate an electric 



field. While EOF and induced-charge EOF is well suited to propel fluid ( Stone et al. 
2004[ [Squires fc Bazant|2004[|Kim et a/.|2002}|Yossifon et al.\2066^ or to enhance mixing 



in micro-channels (Wang et al. 2006 Barz et al. 2011), ion advection is described by 



the so-called streaming potential (Dukhin 1993). It is relevant in a number of physical 
phenomena related to the advection of charged interfaces such as in the electro-viscous 
drag enhancement observed in particle suspensions ( Sherwood|| 1 980| ) . In addition, it can 
be used to convert mechanical (and as will be shown also thermal) into electric energy 
( Yang et aL|[2003) . This is, among others, in the focus of the current study. 

Electrokinetic phenomena are commonly associated with the excess of one ion species 
in the vicinity of an interfacial charge of opposite polarity carried by submerged solid 
bodies or walls. The ions form a diffusion-dominated electric double layer (EDL) which 
screens the surface charge. Unlike the ions in the Stern layer, the ions in the EDL remain 
mobile and, as in electrokinetic streaming applications, can be advected with the flow. 
Depending on the bulk ion concentration, the EDL is typically only a few to a couple 
of hundred nm thick so that many studies of electrokinetic phenomena do not resolve 
the EDL but assume an effective slip velocity (Smoluchowski- limit). The liquid outside 
the EDL is irrelevant for the momentum source term driving the flow. By contrast, it 
contributes to the usually undesired ion flux by means of electro-migration (i.e. the bulk 
conduction current) caused by the existing potential difference. Therefore, to minimize 
the detrimental influence of the bulk fluid, many studies on electrokinetics are concerned 
with strategies to obtain (characteristic) system dimensions of the same order as the EDL 
thickness either by miniaturization (van der Heyden et a/.||2005| |Daguji||2009| |Xie et al. 



2011) or the utilization of (two-phase) plug flow (Lac & Sherwood 2009). In turn, the 



downsizing promotes the importance of interfacial stresses such as in studies related to 
molecular and apparent slip at super hydrophobic walls ( Yang fc Kwok|2004 Zhao|20lT ) 
or free surface-guided electrokinetic channel flow ( Lee et aL|2006 ). These approaches aim 
at reducing viscous losses in areas of the fluidic domain where excess ion predominately 
accumulate, i.e. close to submerged walls carrying a surface charge. Most studies of the 
electrokinetic streaming potential are concerned with pressure-driven flow while compar- 
atively few investigations were performed on shear-driven flow (Song & Wang 2004) or 
other sources of fluid propulsion. In latter cases, owing to the superposition principle 
in Stokes-flow, electrokinetic and fluid mechanical effects can be linearly superimposed. 
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As a result, the streaming potential becomes a linear function of the driving pressure 
difference. 

To date, studies on thermal effects in electrokinetic flow are comparably scarce. Never- 
theless, within the general scope of recent energy sustainability efforts, it is of interest to 
investigate thermally driven electrokinetic charge separation ( Grosu fc Bologa][2010) by 
using waste heat (generated for example by the central processing unit (CPU) of a com- 
puter). Most of the thermally induced fluid propulsion (by buoyancy, thermocapillarity 
or evaporation) can be formulated -at least in the Stokes limit- as an effective pressure 
difference, which can be subsequently combined with the conventional electrokinetic the- 
ory to estimate the streaming potential generated by a thermally propelled liquid. In this 
case, a combined study of thermal, fluid mechanical and electrokinetic effects appears 
not to be necessary. This holds as long as other effects induced by a variation of temper- 
ature are negligible. Roughly, four different non-isothermal contributions may enter the 
problem formulation: firstly, most of the bulk properties such as viscosity, diffusivities, 
electric conductivity and permittivity are temperature-dependent. Under the application 



of direct (DC) (Wong & Melcher 1969) or alternating current (AC) (Gonzalez et al. 



2006) potential differences this may lead, for instance, to electro-convection. Secondly, 



dissipative effects occurring in the bulk, such as viscous dissipation and Joule heating 
dZhao fc Liao|[2002{ [Maynes fc ^\%bb|[2504| [Sadeghi fc S"aidi|[20To| , should be included 
in the energy equation as well. Thirdly, the formation of a wall (C-) potential is strongly 
dependent on the dissociation processes of surface groups and ion absorption at the wall 
(Revil et al. 1999). More specifically, the wall potential is determined, at least under 
quasi-equilibrated conditions, by a temperature-dependent equilibrium constant, leading 
in turn to a temperature-dependent (^-potential. This becomes particularly important 
for flow through porous media at elevated temperatures, as, for instance, treated in geo- 
physical research studies (Ishido et al. jl983, ; Reppcrt fc Morgan 2003). Lastly, similar 



to the thermal diffusion of colloidal particles in a non- isothermal liquid (Piazza 2004 



Wiirger 2010) or the Soret-driven thermosolutal convection in porous media saturated 



with a binary mixture ( LakshmiNarayana et al. 20081, also the charge carriers in an 
electrolyte are set into a thermally induced diffusive motion. Under the condition of local 
charge neutrality, steady-state and the absence of any external pressure gradient, this 
leads to a well-known thermoelectric potential in bulk electrolytes ( Guthrie et aZ.||T949 ). 
As briefly summarized in appendix |Aj this Seebeck-type of thermoelectric potential van- 
ishes if the thermal mobilities of each ion species do not differ from each other. The effect 
is typically quantified in terms of a bulk Soret-coefficient, ctt, and can be enhanced by 
ion-selective membrane technology ( Hills et a/.|1957 Tasaka|1986 ) or by the utilization of 
more exotic electrolytes ( Bonetti et aL||20TT ). It is also relevant for the sensation of heat 
felt by humans ( Tyrrell et al. 1954 ) . Virtually all of these studies involving the Soret- 
effect assume local charge neutrality. While this is valid in the bulk, the fluid within the 
EDL is not electro-neutral. This is particularly relevant for electrokinetic flows through 
non-isothermal nano-channels, e.g. employed as electrochemical thermal energy harvester 
dKang et a/T]|2012[ ). 

This present work focuses on the implications and significance of ion diffusion induced 
by a gradient in temperature (thermodiffusion) in symmetric electrolytes occurring in 
non-isothermal electrokinetic charge separation processes in confined geometry. As a 
model system, a pressure-driven slit channel flow of a fully dissociated binary electrolyte, 
subject to a temperature gradient along the channel axis, is chosen. Unlike the classical 
treatments of thermoelectricity in bulk electrolytes, the condition of local charge neu- 
trality is not enforced. In addition, it is assumed that both ion species have the same 
thermophoretic mobilities (essentially excluding the conventional thermoelectric effect) 
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and -for simplicity- also the same (Fickian) diffusion coefficients. A modified Boltzmann 
distribution is derived from the scaled non-isothermal Nernst-Planck equation, taking 
advantage of the disparate ratio between channel height and channel length. This leads 
to an EDL thickness which exponentially increases with local temperature. While the 
uniform growth of the EDL generally increases the streaming potential, the thermal 
gradient along the channel axis gives rise to a corresponding gradient in the EDL poten- 
tial. In turn, a gradient in electro-hydrostatic pressure and an additional axial Maxwell 
stress acting on the liquid emerges. While the former is known to be responsible for the 



thermo-osmotic propulsion of colloidal particles in a thermal gradient (Derjaguin et al. 



1987 Wiirger 2010) as well as for the thermo-osmotic transport across porous mem- 



branes (Dariel & Kedem 1975), the latter appears to have been commonly neglected. 
It is an elect rohydro dynamic effect and develops due to the combination of the rela- 
tively weak thermophoretic ion motion with the strong potential gradients within the 
EDL. The two effects due to an axial gradient in the EDL potential induce an thermo- 
electro-osmotic flux, either diminishing or enhancing the streaming potential. |Sasidar 



Ruckenstein ( 1982| investigated a similar scenario in the anomalous electrolyte osmosis 



through channels subject to an axial concentration gradient. However, so far and to the 
best of our knowledge, the implications of such a flux on the electrokinetic streaming in 
a non-isothermal microchannel have never been explicitly considered. 

In ^j2j an analytical model based on the Debye-Hiickel (Dll)-approximation at low 
^-potential and a full numerical evaluation of the local, non-isothermal streaming po- 
tential (depending on the pressure and temperature difference, EDL thickness scaled to 
the slit height, ^-potential, Soret coefficient, temperature-dependent viscosity and on the 
temperature-dependent dielectric permittivity) is presented. In f|3j the total streaming 
potential is obtained by numerical integration along the channel axis and discussed for 
a realistic range of parameters. Finally, the thermoelectric potential of a confined elec- 
trolyte in a thermal gradient is obtained and discussed in the limit of a vanishing pressure 
difference. 



2. Model equations and perturbative solutions 

In the following, the governing equations are summarized and simplified to obtain 
analytical solutions of the streaming potential caused by a pressure driven flow inside 
a parallel-plate slit-microchannel of length Iq, and with one half of the gap width de- 
noted by /iQ- The electrokinetic properties are commonly a strong function of the non- 
dimensional parameter k/iq, where = ^J eHg^sToo/ (2e^zy2noo) is the Debye- length of 
the electric double layer (EDL) near charged interfaces. The Boltzmann constant is de- 
noted by /cs, and Too is the ambient reference temperature. With eg being the dielectric 
permittivity of vacuum and £r,iiq being the relative permittivity of the liquid, the (in 
general temperature-dependent) permittivity of the liquid is euq — £Q£r,iiq = f{T). The 
elementary charge is e, v is the valence of the symmetric v : v electrolyte, and rioo is the 
reference volumetric ion number concentration far away from the walls. 



2.1. Double-layer potential 

For sufficien tly low electric fiel ds and ion concentrations distinctively below IM (M = 
mol dm~'^) (Levine et al. 1975), ion transport in liquids is commonly described by the 
Nernst-Planck-equations reading 

dtnk = V • (fn.fcVnfc 



(2.1) 
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where the substantial derivative dt = d/dt is denoted hy dt{.) = dt{.) + v ■ 'V {.) with the 
velocity vector v = {u,w). The volumetric number concentration of positive or negative 
ions labeled with k — (+, — ) is denoted by n^, and T is the local absolute temperature. 
The diffusion coefficients due to concentration gradients are Dn,k, while Dx.k are the 
thermophoretic mobility (i.e. thermodiffusion) coefficients of the ions subjected to an 
external temperature difference AT. The ion mobilities under the action of a gradient in 
the EDL potential, Vf/;, are denoted by ujk ~ Dn.k/{kBT). 

Thermal diffusion in multi-component fluids is frequently described in terms of so- 
called heats of transport, (Helfand 1960), or, equivalently, entropies of transport 
~ Qkl'^ (Tasaka 1986). As summarized in appendix [A] these quantities are accu- 
rately defined within the phenomenological theory of non-equilibrium thermodynamics 
(de Groot & Mazur 1984 Fitts 19621 and emerge from cross correlation between heat 
transport due to material fluxes on one hand and matter transport due to thermal gra- 
dients on the other. Related phenomena are commonly termed heat-matter cross effects. 
The problem at hand is a tertiary mixture of electrically neutral solvent along with two 
ion species, so that, according to the theoretical framework mentioned, cross diffusional 
effects due to concentration gradients are involved. The latter can be omitted if elec- 
troneutrality is assumed throughout the fluid domain, simplifying the problem into an 
effective binary mixture ( Haase|1969 ). By contrast, as particularly relevant for nanochan- 
nel flow, electrolytes close to walls carrying a surface charge are not electrically neutral. 
The problem remains tertiary, thus also involving cross-diffusional effects between dif- 
ferent ion species. This is routinely neglected in most studies of electrokinetic streaming 
involving wall-effects, in which the ion transport is nevertheless described by the (isother- 
mal) Nernst-Planck equation. Adding thermophoretic ion diffusion to this equation in 
terms of might give the misleading impression that the resulting equation is fully con- 
sistent with non-equilibrium thermodynamic theory (needed to define the Q^), although 
this would only be the case if cross-diffusional fluxes between different ion species were 
included. To avoid this source of confusion and to emphasize the (still) limited validity of 
(2.1), instead of effective thermal diffusion coefficients Z^T.fe are used herein. Never- 
theless, within the present approximation, the distinction between DT,k and has solely 
a cosmetic character. In fact, in appendix [a] it is shown that (herein) 5^, = Q*^/ {ksT'^) 
( Wurger||2010 1, where = DT.k/Dn.k are the intrinsic Soret coeffic ients (other a uthors 
would call it the thermal diffusion ratio) of the ions in units of ( Vigolo et al. ]|2010| ). 
In this context, it is important to point out that all of these parameters equivalently 
quantifying the thermomobility of individual ion species (Q^, Dx^k or S'^) can be experi- 
mentally determined only relative to each other ( Hills et aL|1957 ) and not on an absolute 
scale. 

In the following, it is assumed that the flow is stationary and fully developed so that 
(it rife = udxUk- To estimate the characteristic axial flow velocity Mq, used in the following 
to scale the governing equations, a conventional parabolic velocity proflle in a slit channel 
is employed so that ug — '-^PaAho/{3rio), where A = ho/h < 1. The external pressure 
difference is denoted by Apa and rjo = ri{Toa) is a reference dynamic viscosity of the 
electrolyte. For sufficiently viscous fluids and sufflciently small values of Apo and Hq, uq ^ 
0{A), i.e. uo can be expected to have low values. This assumption will be verified later 
on. Therefore, the common assumption is made that the problem is diffusion dominated. 
In fact, if the ionic Peclet-numbers Pen.k — loUo/Dn^k are 0(1), then the ratios between 
the characteristic times of ion diffusion, t^^k — h^/Dn^k, and the characteristic advection 
time, tu = Iq/uq, are 0{A^). In this case, the substantial derivative in (2.1 1 can be 
neglected to first order in A. In addition, while temperature-dependent particle mobilities 
are well verified in isothermal flow, under non-isothermal conditions, they essentially 
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violate the linear Onsager reciprocal conditions. The present analysis is based on the 
general framework of linear response theory of thermodynamic non-equilibrium processes. 



Thus, in (2.1), the particle mobilities are assumed to be not depending on any of the 
thermodynamic conjugate driving forces, including temperature. Particularly, herein uj = 
Dn.k/ (kBToo)- It follows that 

VNk + NkSlATVe + VkNkV^ « (2.2) 

with Nk — nk/rioo, Q ~ (T — Too) /AT and — evtp / {kBToo)- The -D„,fc must not be 
zero, Vk = Vkjv and c\ is an integration constant. The latter is set to zero, implying that 
all fluxes vanish if no gradient in number concentration, temperature or EDL potential 
is present. Consequently, the local ion number density N-^ follows the apparent electric 
EDL potential — + S^ATQ according to the Boltzmann distribution Nk = 
C2exp(— ^ffe). Thus, in case of non-isothermal conditions, the electric potential needs 
simply to be corrected for the intrinsic Soret effect. The integration constant C2 is set to 
unity, i.e. all ion species have a uniform density rioo at = 0. This corresponds to a 
uniform thermo-electro-chemical potential, /i^ = v^eip + kBToo[^^ink) + S^T], throughout 
the domain. According to this derivation, the relationship should be valid for diffusion 
dominated conditions and as long as the 5*^ are constants, while Dn,k may vary e.g. with 
temperature. For a symmetric v : z/-electrolyte with K species, the charge density pf 
reads 

K 

= y^VkNk = -exp(-S'* ATe)[exp(*) - exp(-AS'*Are)exp(-*)], (2.3) 

°° k=l 

where AS** ~ S'/^ — S't. For small differences in the intrinsic Soret-coefhcients of the 
anions and cations, S"^ w 5*1 = S* and small temperature differences AT, it follows 
exp(— AS"* ATO) « 1. In this case, the Soret effect only diminishes the local charge 
density by a factor of /|. = exp(-S'* ATO). Then 

K 

-^-^]7feiVfe = -2/|.sinh(*). (2.4) 
eiyfioo , ^ 

k — l 

Note that, in the present case, AS"* = infers identical mobilities for both ion species 
when exposed to a temperature gradient. Therefore, it has to be emphasized that the ef- 
fect conventionally causing a thermoelectric potential is explicitly excluded in this work: 
For instance, the Seebeck effect observed in electric conductors appears due to different 
mobilities of the charge carriers. This is especially apparent in semi-conducting solids, 
where the electrons populating the conduction band have a much higher mobility than 
the positively charged donors left behind during the temperature-dependent population 
step. The difference in the charge carrier mobility leads to the observed electric poten- 
tial upon exposing the material to a temperature gradient. On the contrary, here it is 
assumed that both ion mobilities in a non-isothermal matrix are the same, i.e. no ther- 
moelectric potential should be generated. Nevertheless, in a later section it is shown that 
a thermoelectric potential is generated even under the present condition as a confinement 
effect. 

The charge density is related to the EDL potential as expressed by the Poisson equa- 
tion, namely V • (eVV') = —Pf- The axial x-direction is scaled with Iq while the lateral 
z-coordinate is scaled with ho-, i.e. X = {X,Z) = (x/lo,z/ho)- It is assumed that ho is 
much smaller than the axial extent lo and thus <C 1. In non-dimensional form, the 



Flow and streaming potential axial, non-isothermal 



7 



Poisson equation reads 

+ E*ATdx-^dxe) + + E* /STdz-^dzQ = - ^I7fciVfe, (2.5) 



K 



k=l 



with 7t = k/iq and d^p = d/dtp (tp = X,Z). The temperature dependence of the dielectric 
permittivity was incorporated with Ve = et'VT where ct = dcuq/dT and E* = CT/euq- 
The terms with derivatives in axial direction are small to order and can be neglected. 
The temperature distribution is governed by the energy equation. Neglecting viscous 
dissipation, Joule heating as well as the kinetic energy of the flow and assuming a constant 
thermal conductivity fchg.o of the electrolyte, this is approximated by dtT = agV^T. The 
thermal diffusivity is denoted by = q/{cpPq), where Cp is the heat capacity at 
constant pressure and pg a constant reference fluid density. In non-dimensional form, the 
energy equation reads 

APe {drO) - A'^dj^e = 9|e (2.6) 

where Pe = hQUo/ao is the thermal Peclet number. In non-dimensional notation, the sub- 
stantial derivative is dr{-) = dr{-)+V-'V {.) = {lo/uo)dt{-), with t = tuo/lo, V = {dx,dz) 
and V — {U,W) — {u/uo,w/wo). As implied by the non-dimensional continuity equa- 
tion, V • V, the vertical scaling velocity is wq = Auq. For electrolytes with Pe ^ 0{A), 



the left hand side (LHS) of (2.6) can be neglected and -using the symmetry condition 
along the centerline a.t Z ~ 0- it follows that the temperature is identical to the local 
wall temperature. It is assumed here that the wall temperature varies linearly in X, so 
that also the liquid temperature varies according to dxQ — constant, while dzQ = 0. 



With relations (2.4) and (2.5), this leads to 

al^f = K|.sinh(*) (2.7) 

where ks* — ~i^fs* = KhQexp(—S*ATQ/2). It follows that, within the limit of this ap- 
proximation, the Soret-effect simply implies a local, temperature-dependent thickness of 
the EDL according to k^} = (k/s*)"^; on the other hand, the EDL-potential distribu- 
tion is qualitatively the same as for isothermal conditions. The case of dzQ = constant, 
while dx& = was treated in a separate publication ( Dietzel fc Hardt|2012 ) 



For comparison, if one disregards the Soret-term in the ion transport equation (2.2) 
{S'l — 0) while keeping the temperature-dependence of the electromigration mobility, one 
can deduce that the EDL thickness is modified according to = k^^ ^/T+~AT/T^. 
As an order of magnitude approximation, it is assumed that S* ~ I /Too • This leads to 

_ n,} ^ exp(e/2) 



where ^ — AT /Too- A Taylor-expansion of (2.8) gives x ~ 1 + C /4- This indicates 
that, for systems violating the linear response theory (unlike those treated in this work), 
the thermophoretic (Soret)-term in the ion transport equation and the temperature- 
dependence of the electro-migration mobility may expand the EDL-thickness to a similar 
extent, especially for small AT and small S* . 



The solution of (2.7) has been extensively discussed in the literature. Within the com- 
mon Debye-Hiickel (Dll)-approximation [^E* <C 1 — >■ sinh(\I') vp], the EDL-potential 
distribution reads 

^nn^C^s'^^^^^P^- (2.9) 

' COSh(K5. ) 

The dimensionless ^-potential at the wall is denoted with C,\g — evC\s/{kBToo), where C,\s 



8 



M. Dietzel and S. Hardt 



is the original ^-potential. For completeness and later reference, the main steps leading 



to the full solution of (2.7) and comparison with (2.9) are repeated in appendix B 



As it becomes relevant and will be discussed in the next section, the axial gradient of 
^ does not vanish. In this work, a constant ^-potential along the channel is assumed. 
Ignoring the charges in the Stern-layer and denoting the axial gradient of the EDL at 
the wall with dzip\s7 the surface charge q'/ = —eugdzipis — —QiqkBToo / {hQeiy)dz'^ \s can 



be calculated with (2.7) to read 



q\s = fs' \/S>QiqnookBToo 



\ 



sinh 




sinh 



2 I *|c 



(2.10) 



Here, the (non-dimensional) EDL-potential at Z = along the channel axis is denoted 
with — evip^cl {kBTao, where is the original center-line EDL-potential. Under the 



condition of non-overlapping EDLs (5'|c = 0) and isothermal conditions [fs* — 1), (2.10) 
is identical to the classical Graham-equation. Herein, "^^^ and /g. are both functions of the 
axial coordinate. Consequently, at constant Cj^, q'^^ is not constant along X. Alternatively, 
if one assumes a constant surface charge along the channel walls, the EDL-potential 
according to the DH-approximation is given by 



DH.q'' 



q'\s cos\\{ks'Z) 
fs* sinh(K5-) 



(2.11) 



with g"^ ~ q'^s/ ^/'^eiiqUookBToo- Expression (2.1l|) implies that dx^ still does not vanish. 
Thus, the following derivation is quite general. 



2.2. Axial velocity profile 

The motion of an incompressible electrolyte of density p = pq and mass averaged velocity 
vector V is described by the Navier-Stokes equation with Maxwell stresses added as a 
source term, namely pdtv = V • cr^ -I- V • tr*^. The mechanical stress tensor is cr^ = 
~pl + ril'Vv + (Vw)"^] with / being the unit tensor, and where the dynamic viscosity rj 
does not necessarily need to be constant but might vary with temperature. The small 
channel size allows to omit hydrostatic contributions. For a homogeneous, incompressible 
fluid, the Maxwell stress tensor reads ct^^ = eiiq{Vcf)V<j) — Vcf)-V(j)l /2). According to the 
Osterle-postulate, the total electric potential, = V + 0sti is a linear superposition of the 
electric potential due to the ion double layer at the interfaces, V', fulfilling the Poisson 
equation, and an externally applied electric potential, (j)st, complying with 'S/'^4>st = 0, 
i.e. being source-free within the electrolyte. The force contribution due to the Maxwell 
stress tensor can be expressed in terms of the Kortew eg-Helmholtz electric fo rce per 



PfV(j)- (V</>)2eTVT/2 (Russell k Saville 1989). With 



volume, namely V 

$ = eiy(j)/{kBTao), the non-dimensional axial velocity component fulfills 

ARe {drU) - A^ {r]d\U + 2dxUdxf] + dxWdzrf) 



-A^^{dl<l>dx^ + E*/S.T 



dx^dx<S> - \ (9x$)' 



dxQ \ = -dxP 



-dz ivdzU) + 1^ <; $ + e*at 



dz-i'dx<i>dze^^{dz<^>fdxe 



(2.12) 



while the lateral component is obtained from 

A^Re (drW) - A^ {A^rjO^W + rjdlW + A^dxWdxV + 2dzWdzr] + dzUdxv) 
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-A^^\dlmz<f + E*AT 



+ 5 <! d'z-^dz^ + E*AT 



dz-^dz^ - \ [dz^f 



= -dzP 
dzQ \ ■ (2.13) 



In this formulation, the charge density pf was expressed with ( |2.5[ ), Re — puoho/rjo is the 
Reynolds-number with rjQ, and rj = rj/rjo is the non-dimensional local viscosity. Further, 
Ha = 2AhQnaokBTao / {uorjo) is the (scaled) Hartmann-number, and the fluid pressure 
is non-dimensionalized according to P = Ah^p / [uqI^q) , with p being the original pres- 
sure. As discussed earlier, the characteristic velocity uq is a linear function of A and it 
was assumed that wq = 0{A). This suggests that the i?e-number can be expected to be 
small, being at least of order A as well. In addition it follows that Ha (as well as P) 
does not depend on A. Since the ionic Peclet-number Pen,k was already assumed to be of 
0(1), Ha cannot be of the same order ( Yariv et aL||201l" |. In fact, also with the present 
scaling one has Ha = (<j/Peji_fc)K^, where ^ = eiiq/(riQDn)[kBTaol{ve)Y' is the intrin- 



sic Peclet-number (Saville 1977[ ). For typical aqueous solutions ~ 0.5. Consequently, 
Ha/K^ ^ 0{1) is a consistent, and for the present purpose sufficient scaling. All of these 
assumptions need to be verified later on. Furthermore, it is assumed that the external 
field 4>st is caused by convective ion streaming in axial direction and only dependent on 
the axial coordinate so that dz^ = S^Vf (and thus, = 9^5' as well). Also, the ther- 
mal gradient has only a n axial component. Neglecting terms of order A^ and smaller, one 
can deduce from (2.131 that P k, i/a/(2K^)(9z^')^ -I- c^{X) where c^{X) is an integra- 



tion constant. Here, the latter is simply the externally applied pressure Pa{X). The first 
term proportional to the square of the lateral electric field is the electrostatic pressure 
contribution, representing the electro-osmotic pressure of the ion cloud, Posm = nkBT, 
n — X^s^i ^fc- common studies of electrokinetic streaming in long microchannels, this 
term is not a function of the axial coordinate direction X. By contrast, as shown before, 
the thickness of the EDL and the EDL-potential are herein a function of temperature due 
to the Soret-effect and the change of the electric permittivity with temperature. There- 
fore, the electrostatic pressure varies in axial direction, and the overall axial pressure 



gradient is dxP = Ha/ {2K^)dx{dzW + dxPa- With Re 
the LHS of (2.121 is small to order A^ and one can write 



0{A) and Ha/n^ < 0(1), 



dz ivdzU) « ^ 



dx (dzW ~ dz (9z*9x$) 



]^{dz-^f E*ATdxQ 



^dxPo- (2.14) 



In general, the viscosity is a function of the shear rate, concentration of dissolved species 
as well as of temperature. Shear rates are assumed to be sufficiently small so that shear 
thinning or thickening behavior is of no importance. Furthermore, significant relative 
concentration changes of the dissolved ions are only present in the EDL. Absolute val- 
ues of ion concentration are proportional to rtoo which is typically very small in dilute 
electrolytes as treated herein. Therefore, even within the EDL, the dependence of the 
viscosity on the local ion concentration is expected to be negligibly small. Finally, the 



temperature varies only in axial direction and (2.141 can be integrated twice in Z even 



without explicit knowledge of the viscosity-temperature-correlation. Symmetry holds at 
Z — (subscript c), while the no-slip condition has to be fulfilled at the wall {Z = 1, 
subscript s). With (f> = 4* + ^st, this leads to the expression for the axial velocity, where 
= ev(l)st/{kBToo)- 



U ■ 



dxPg 
2f] 



{l-Z^ 



Ha 



9x<fst(*-C|J 



Ha 

9 



[n~Vt\s + dz^\c{l- Z)]. (2.15) 
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Note that the local dimensionless viscosity might still depend on the axial coordinate, 
i.e. fj = f(X). In expression (2.15), one has 



dx {dz'^f d^Z - / dz'fdx'fdZ 



dxQ 



E*AT I I {dz'i'fd^Z, (2.16) 



where the integration symbols denote primitives of the corresponding functions and d^Z 



denotes double integration of the integrand with respect to Z. In (2.161, the first two 



integrals are the sum of the electrohydrostatic (EHS) contribution and the electromigra- 
tion force (EMF), while the EHS alone is just one half of the first integral. For symmetry 
reason, it is 9z^|c = and thus also dz^\c — 0. The expression for fl vanishes if the 
EDL-potential 5" does not depend on the axial coordinate and a temperature-dependence 



of the permittivity is neglected. In this case, equation (2.15) resembles the well-known 



isothermal result. In this work, depends on the local value of the Debye-parameter 
^S' = f{^)- Thus, dx'^ = chist'^dxKs' , where it follows from the expression for Kg* 
that 

dx{KS') 1 



Kg. 



-ATdxO{S* +E*). 



As derived from (2.7), dz"^ = ±2 Tcs- [cosh^(^'/2) - cosh^(^'|c/2)]i/2 and thus 

dxKs- 



dx {dz'^f = 2 {{dz^f + 4. [sinh(*)as,. * - sinh(*|,)ft.,. ^-i,] } 



US' 



(2.17) 



(2.18) 



The derivative 9k„. '5 can be evaluated from the implicit equation (B 2 ) e.g. by bisection 



or an iterative fix point procedure. Yet, even for small values of k^*, for which does 
not approach zero, b is nevertheless close to unity. For this value of b, F(7r/2,fe) has a 
singularity, leading in turn to difficulties in accurately evaluating the implicit function. 
Therefore, in practice, it is preferable to evaluate this gradient numerically. For this 
purpose, (2.7) is solved for several values K5. — ks',w in the vicinity of a mean value 
^S*,o, where the difference ks',w — '*s*,o = Ak is a multiple w of a small increment 
Ak <C 1. For instance, the derivative dji^, ^' can be approximated with a finite difference 
(FD)-scheme, using the discrete distributions of 4'(Z)|jjg, ^ solved in the previous step. 
Alternatively, assuming an EDL-potential based on the DH-approximation as derived 
before, one has 

i^s' '^)dh = *CH [tanh(K5. Z)Z- tanh(K5. )] . (2.19) 
The partial derivative ^^^'(Z)/^!^ = 9k-0(Z)/C|s is plotted in figure [T] for several values 
of the dimensionless Debye-parameter k = (0.5,1,5,20). In this context, to indicate 
universal applicability, ks" was substituted by k, and the plots are valid for any functional 
correlation k{X). The wall ^-potential is either C|s = 25 • 10^'^V, corresponding to C|s ~ 1 
and shown in (a), or C,^g — 125 • lO^'^V, corresponding to ~ 5 and shown in (b). 
The dashed (red) lines marked with a symbol were numerically evaluated, employing 
a 5-point stencil FD-scheme, which is accurate to third order in Ak (Fletcher 1991[). 
After a grid independence study, Ak = 10"'^ was chosen. As detailed in appendix |B[ 
(2.7) was solved for each Kg- ^ by collocation with the BVP4C-function implemented in 



Matlab. Approximately 1200 grid points in Z-direction were used to obtain solutions with 
a relative tolerance of 10~*. For comparison, solutions based on the DH-approximations 
were computed with (2.19) and are shown in Figure [l] with (green) dotted lines. In all 
plots, the different line symbols indicate the value of k used. At the wall [Z = 1), one 
always has 9^^' = since a constant C-potential is assumed. For both values of C,\s and for 
7c ^ 1, the local minimum of i9k5'/C|j, is along the center line {Zmin — 0) while for larger 
values of Tc, the minimum is located in the vicinity of the wall, with Z^in > 0-8. Absolute 
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Figure 1. Partial derivative ftf*(Z)/C|^ = ftr?/'(Z)/C|3 with k = (0.5,1,5,10). (a) 
C|s = 25 • 10~^V (C|s ~ 1), (b) C|s = 125 • 10"^V (C|, ~ 5). Thick dashed (red) lines mark the 
numerical solution (NM), while results marked with thin dashed (green) lines were obtained, 
within the Debye-Hiickel-approximation (DH), from (2.191. 



values of 9k5'/C|s decrease with increasing Debye-parameter and practically vanish along 
the channel center line for k ^ 5. In these cases, ^'|c ~ 0, leading to a vanishing derivative 
Sk^E' as well. For large Tt, the DH-approximation captures reasonably well, even 

for — 125 • 10~'^V. On the contrary, for overlapping EDLs, the DH solution locally 
deviates from the numerical solution by up to 50%. For C|s = 125 • 10~^V and k ^ 1, the 
numerical solution is steeper close to the wall and flatter along the center line than the 
DH-solution. 

After some algebra and within the Debye-Hiickel-approximation, one finds for the 51- 



integral expressed by (2.16) 



-cosh(2K5*Z) 



tt2 72 



8cosh^(K5* 

Since dzi^\c = 0, the axial velocity distribution becomes 

Udh = - Z ) - __2 Ox'^st 



S* + 2Kl,Zhanh(jis')iS* + E* 



(2.20) 



2?7 



rj K 



cosh(Ks-^) 
cosh(Ks* ) 



- 1 



Ha C\s^Tdx& 
rjK^ 8cosh^(Ks*) 



cosh(2Ks.Z) - cosh(2K5.) 2 n r^2^ 



S* 



2 

27t|.(l - Z2)tanh(Ks.) {S* + E*)] 



(2.21) 



Expression (2.21) describes the axial velocity profile across the channel if, next to ex- 



ternally applied gradients in axial pressure and electric potential, a thermal gradient 
is present along the channel as well. It is applicable to any type of electrokinetic flow 
(EOF or SP) in a slit-micochannel. It is emphasized that U (and Udh) denotes only 
the convective fluid motion. As usual, any (axial) diffusive flux induced by temperature- 
or concentration gradients is relative to U (respectively Udh)- Nonetheless, under the 
present assumptions and non- isothermal conditions {dxQ 7^ 0), an axial convective flow 
is -counterintuitively- present even without a pressure difference {dxPo = 0) and without 
an external field {dx^st — 0). Under non- isothermal conditions, the combination of the 
-weakly- varying EDL-thickness along the axial coordinate X due to the Soret effect (and 
^liq = f{T)) with the steep lateral gradient of tp inside the EDL leads to an axial field in- 
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side the EDL. This field causes an axial gradient in electrohydrostatic pressure as well as 
an axial electro-migration force, which both set the fluid in axial motion. Hence, the non- 
vanishing axial velocity under non-isothermal conditions is essentially an electro-osmotic 
flux. Under the common assumptions made in conventional isothermal treatments, one 
can prove (see appendix [C]) that -as expected- the EDL-field itself cannot propel the sur- 
rounding fluid since the EDL-field is always perpendicular to the streamlines. Unlike for 
isothermal conditions, an axial temperature gradient removes this orthogonality and the 
EDL-field drives the flow. As shown in appendix [C] this is not an artifact of the lubrica- 
tion approximation or the specific scaling used in the present derivation. As discussed in 
a later section, this thermo-electro-osmotic flux also induces a streaming potential. This 
does not contradict the second law of thermodynamics, as a constant heat flux through 
the channel is required to maintain the temperature gradient. It is also to be empha- 
sized that, in this study, the incorporation of the Korteweg-Helmholtz electric force in 
the momentum equation is of crucial importance; otherwise the described phenomenon 
cannot be correctly captured. This distinguishes this study from common isothermal 
electrokinetic considerations where the force term is not necessary to predict a streaming 
potential but is merely included to consider the (commonly weak) electro-osmotic flux, 
opposing the pressure-induced flow, and to fulfill the Onsager reciprocal condition. 



2.3. Electric Currents and Streaming Potential 

At stationary conditions, the total electric current due to ion motion driven by advection, 
concentration- and temperature-induced diffusion as well as electro-migration vanishes. 
From this definition, the streaming electric field E^t — —dx^st in axial direction can 
be derived. The total axial streaming current is 1st — 2 Ay JI^° pfudz with Ay as the 
extension of the micro-channel in y-direction. This current is not necessarily constant 



along X, i.e. in general 1st = lst{X). According to the leading order terms in (2.5) the 
charge density pf/{evnoo) = J2k=i'^kNk is equal to — 2/(k^)9|5'. Integrating by parts 
under the given boundary conditions leads to 1st — '^■AyuQekBTrx,/ QiQev) dzUdz'i'dZ . 
With an expression for dzU derived from (2.15), this leads to 



Istve 



2AyKUoekBTo, 



where 



st.Pn 



1 



Ist.PodxPa + Ist.K.t^x^st + Ist.edxQ, 



■^dz] , Ist.^^, = / {dz-^fdZ, 

/ 1^ V Ja 



(2.22) 



Ist,e — 



Ha 1 



K rj Ox 



dxO 



dz^dz-^dZ. 



(2.23) 







With (2.16) and (2.18), the last integral in (2.23) can be written as 

r-l 



dz^dz'^dZ = -Ks- 







dxus* 



US' 



[idz^fdj,^,^}dZ 



dxKs- 

Ks- 



dz^ / [dzW 



1%. [sinh(*)as,. * - sinh(*,)9^^. ■^,] dZ j dZ 
E*ATdxe , I ^g^^y^z^dZ, (2.24) 



where dxKs* / ks" is given by (2.17). In general, "^{Z) and (^^.^'(Z) are obtained from 



numerically solving (2.7) (see previous section). The derivative of the EDL potential in 
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Z-direction is given by dz'^ = 2 7is*[cosh^(*/2) - cosh^(*|c/2)]i/2. Thus, the integrals 
in (2.23) and (2.24) are fully described and can be numerically evaluated. 

With expression (2.21), = '*s*C|sSi'^h(K5«^)/cosh(K5*), and after some algebraic 
manipulations, the (dimensionless) streaming current can be written within the DH- 
approxiniation as 



2AyKUQekBToc 
2k r] 



tanh(K5'.) 



HaCy,dx^st ^ 
2k f] 



tanh(K5* 
Ks- 



cosh {ks') 



tanh(K5*) tanh'^(K5*) Ks.tanh(K5. ) — tanh^(K5. ) — 1/2 



27t5. 



3ks' 



cosh^(K5. ) 



S* 



E* } (2.25) 



KS»tanh(K5. ) — tanh (ks*) — 1/2^ 1 

cosh^(K5.) 2cosh^(Ks.) 

In steady-state, the streaming current is balanced by currents induced by axial gradients 
in ion concentration, temperature and electric potential. For the case that the coefficients 
of thermodiffusion, respectively of Fickian diffusion, are the same for each ion species, 
the total current resulting from these effects, for simplicity termed "diffusion current", 
reads 



'2AyD„ 



dxPfdz + S*dxT / pfdz 



B-Loa JQ 



K 



dx4>y^^vln],dz 



k=l 



(2.26) 

Similar to 1st, in general Idf = Idf{X). As mentioned, the charge density can be e xpre ssed 
with pf/{evnoo) = -2/{K^)dl'^. Furthermore, Oxk/k = -E*ATdxO/2. With ([zt]) and 
a number of algebraic manipulations one finds 



hf = 4.AyiyenooAD„ [k-^ [dx (9z*)|, + E*ATdxe (^z^-; 



+K-''S*ATdxQ {dz^s - f's'dx [Tig'' (^z*),,] - fs'dx<^st / cosh(^)dz| (2.27) 
One has 



fs^dx 7«^?(9z*)|, ^K-^dx{dz^)^s-2K-^dz^s ^ 



dx (ks*) 



KS' 



(2.28) 



yielding 



Idf ^ ^Ayven^ADr, \ k ^ [dz^) 



[S* +E*)ATdxQ + 2 



dx {ns') 



KS' 



-fs'dx<^stj^ cosh(^')dz|. (2.29) 



According to (2.17), the term in the square brackets is identically zero. Hence, the re- 
maining diffusion current resembles the conventional conduction current due to electro- 
migration induced by the external electric field, namely 



Idfve 



with 



2AyKUQekBTo, 

Dr. 



UqIq' 



fs'Ks' / cosh(*)dZ. 



(2.30) 



(2.31) 
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It is noteworthy that while the streaming current directly depends on terms proportional 
to E* as well as S*, the corresponding terms in the diffusion current exactly cancel. In 



the present study, the first three terms on the right hand side (RHS) of (2.27) are exactly 



compensated by corresponding terms stemming from the electro-migration due to the 



axial EDL-field (fourth term in (2.27)). Based on the given assumptions, this result is 
expected: As derived in section |2.1[ the thermo-electro-chemical potential of each ion 
species equals /i^ = vj^eip + kBTao[l'n{nk) + S* kT] and is uniform throughout the channel, 
while diffusive fluxes can only be caused by gradients in chemical potential. Consequently, 
the former vanish in this study with the exception of the one caused by the external field 
since the latter is assumed to have no effect on the local charge density. Hence, in this 
study, thcrmophoretic ion motion as well as a temperature-dependent electric permit- 
tivity influence the diffusive current only indirectly through the temperature-dependent 
The DH-approximation yields: 



EDL-thickness k^} 



ihf) 



DH ' 



2A?/KUoefc_BTo< 



Dn 
UqIq 



(is 



tanh(K5* ) 



1 



KS' 



cosh'^(K5.) 



(2.32) 



The streaming potential is calculated by setting the total current, being the sum of the 



streaming current [expression (2.22)] and the total diffusion current [expression (2.30)], 
equal to zero. After reinserting the definitions of the dimensionless parameters to obtain 
dimensional values one finds for the local streaming field per pressure gradient 



'x<fst,p 



ksToo An ht,p„+IstMdxO/dxPo 



(2.33) 



Note that while the characteristic parameters A and uq need to stay in a certain range 



so that the simplified governing equations remain valid, the value of equation (2.33) does 



not depend on the specific choice of these parameters. This is because the dimensionless 
parameters Ist,ip = Po, ^st, ©) and /d/,*^^ depend on these parameters in a reciprocal 
fashion in comparison with the pre-factor so that the specific values cancel out. This 
can best be observed by considering the local streaming potential in the DH- limit, which 
reads 



'st,p 



dxPO / DH 

tanh(K5» ) 



Ks- 



1 ^liqC\s 



1 ^iiqC\s 
1 

2cosh^(K5.) 



tanh(Ks'. ) tanh'^(K5* 



27ts. 



E*AT 



Fes 



3ks 

1 <^liqC^! 



S*AT 



2 DnT] 



Fr, 



(2.34) 



where C|s is the dimensional surface potential, AT/Apo = dxT/dxPo (the temperature 
and the pressure gradients are constant herein), 



Fes = I cosH^)dZ « 1 



<\s 



tanh(K5.) 



and 



tanh(K5.) 



cosh (ks* 



cosh (ks* 

Ks*tanh(Ks*) ~ tanh^(Ks*) ^1/2 
cosh^ (ks* ) 



(2.35) 



(2.36) 



As mentioned before, equation (2.34) indeed does not depend on A or uq. It is a first order 
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approximation in A and requires that < 1, AT/T^ < 1, Pen,k ^ 0(1), S'^AT ^ 
0(1), Re < 0{A), Pe ^ 0{A) and Ha/n'^ 0(1). It remains valid if the local values of 
the viscosity and of the diffusion coefficient vary in axial direction, e.g. due to a variation 
of temperature. To limit the scope of this study, this effect will be neglected in the 
following while the temperature-dependence of euq is kept. In the particular case treated 
herein, the temperature-dependence of rj and of £)„ has only a parametric influence on the 
streaming potential. On the contrary, a temperature-dependent dielectric permittivity 
involves additional terms in the Poisson- and momentum equations, leading to effects 
comparable and as interesting as observed for the intrinsic Soret effect. 

For a vanishing temperature-difference, (2.341 agrees with the well-known expression 
of the streaming potential for slit-micro-channels with small (^-potentials discussed in 
the literature (Masliyah & Bhattacharjee 2006). For a nonzero temperature difference 
though, the streaming field according to (2.33), respectively ( 2.34[ ), is not a constant but 
a function of the axial coordinate x. The total streaming potential difference per applied 
pressure difference along the length of the channel, A(j)st,p/ Apo, has to be computed 
by numerical integration of (2.33), respectively of (2.34|, from ^ x ^ Iq. In the next 
section, the expressions (2.331 and (2.34 1 are analyzed based on realistic thermophysical 
properties of the electrolyte and under conditions fulfilling the assumptions made in the 
course of the derivation. 



3. Analysis of specific cases 

3.1. Streaming potential 

To qualitatively assess the contribution to the charge separation of one effect (pressure- 
induced streaming 'P', Soret-effect 'S', temperature-dependent permittivity 'E') relative 
to another, the following dimensionless numbers are defined, which are ratios between 



the leading pre-factors of each contribution appearing on the RHS of Eq. (2.34 1: 



Defining these ratios, the largest possible absolute magnitude of each contribution is 



evaluated since the functions in Eq. (2.34) depending on ks* are 0(1) or smaller. Com- 
mon values for these ratios are listed in the last two columns of Table [T] In this table, 
the electrolyte properties are summarized in the first two columns, where the solvent 
properties are based on pure water and given in the first three rows. The electric and 
transport properties of the solute, listed in the fourth to seventh row of the first two 
columns, refer to a 0.01 M NaCl-electrolyte solution. The corresponding values for a 
KCl-electrolyte are of the same order of magnitude. If not stated otherwise, all val- 
ues were determined at 25°C (Too — 298K). In this study, AT, Apo, Cs, '^^ — '^^o 
and S*AT were varied. The ranges of values taken by these parameters are summa- 
rized in the third and fourth column. The fifth and sixth column provide a selection 
of calculated parameters, which are relevant for the verification of the scaling used in 



the derivation. In this context, A — 0.1 and Uoo — 0-OINa (Mansouri et al. 2007) 
{Na is the Avogadro number) were fixed values so that the reference EDL thickness 
is calculated to be « 10~''m. The nominal channel height is determined accord- 
ing to ho = k/k, while the channel length is determined with Iq — k/(Ak). As men- 
tioned, a conventional parabolic velocity profile in a slit channel was assumed to esti- 
mate the characteristic flow velocity, i.e. uq = ApoAK/{3Kr]Q). The Reynolds-number 
is then Re = ApApQK^ / {Sk^tJq) , the Peclet-number Pe = AApqk'^ /{3K^aor]o) and the 
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Fluid properties || Varied parameters || Scaling parameters || Char, ratios 



po(kgm ^) 


9971^1 


AT(K) 


5-25 


ito(ms ^) 


10^^-10"^ 


S/P (-) 


10^^-10"^ 


?7o(mPa s) 
ao(m^s"^) 




|Apo|(Pa) 


lO^-lO^ 


Pe (-) 


10-^-10° 


E/P (-) 


10-3-10^ 


1.45 ■ 10~^[^1 


C|s(mV) 


5-125 


Re (-) 


10-^-10"^ 


S/E (-) 


10^^-10^^ 




10-9 [3] 




10-^-10^ 


Pe„.k (-) 


10-^-10^ 










S'AT (-) 


10~^-10"^ 


Ha/K^ (-) 


10-4_io2 






^Uq/SO (-) 


78.14^1 





























Table 1. Thermophysical properties and characteristic numbers. The electrolyte properties are 
listed in the first two columns, where the solvent properties are based on pure water and listed 
in the first three rows. The electric and transport properties of the solute, listed in the fourth 
to seventh row, refer to a 0.01 M NaCl-electrolyte solution. The corresponding values for a 
KCl-electrolyte are of the same order of magnitude. If not stated otherwise, all values were 
determined at 25°C (Too = 298K). The third and fourth column list the range taken by the 
parameters varied in this study. The fifth and sixth column provide a selection of calculated 
parameters, which are relevant for the verification of the scaling used in the derivation. In this 
context, A = 0.1 and rioo = O.OINa are fixed values, leading to a reference EDL thickness of 
^ lO^^m. The nominal channel height is given by ho = k/k, while its length is determined 
with lo = K/{Aii ). The seventh and eight columns provide the characteristic ratios expressed in 
|3T|. [1] -[Lidel([200 9), [2] - Buclmcr et al. (1999), [3] - Takcyama & Nakashima (1983), [4] - 
Takeyama fc Nakashima| ( |1983p , iAgar fc Turner| ( |1960p , |Leaist| ( il990| 



Hartmann-number Ha = GnoofceToo/Apo- The ionic Peclet-number, Pen,k = lo^o/ Dn^k, 
as the ratio between the advective and diffusive transport (based on the corresponding 
time scales) is Pcn.k = ApoK^/(3K^-Dn?7o)- As summarized in the table, the selected 
pressure differences are in the range of IPa ^ ApQ ^ lO^Pa, which corresponds to a pres- 
sure gradient of lO'^Pam"^ ^ l^j^pol ^ lO^Pam^^. For comparison, pressure gradients 
typically appUed in studies concerned with electrokinetic streaming are in the order of 
|a^po| = 0( 106) -C>(109)Pam-i ( |Yang aL|2003||Mansouri et aL|2007l|van der Heyden 



\et al. 2005). As mentioned in section |2j the intrinsic Soret-coefRcients (or equivalently, 
the ionic heats of transport of each ion species) can be experimentally determined only 
relative to a reference ion ( Agar fc Turner||1960 ), where the thermophoretic mobility of 
the latter is arbitrarily set to zero. Herein, it is specifically assumed that all the S^. (re- 
spectively the Ql) are identical for each ion species k, i.e. one cannot specify a reference 
ion. Hence, one may rightfully argue that it is in fact not readily possible to determine S* 
experimentally. For this work, this circumstance is ignored and the order-of-magnitude 
of these coefficients provided in the literature are taken as a rough estimate. 

From Table [T] it becomes apparent that the contributions due to the Soret-effect to the 
electrokinetic streaming are small compared to the pressure-induced contribution even for 
these relatively small pressure differences Apo- Therefore, the contribution is insignificant 
for large Apo, as it is typically the case for electrokinetic streaming applications. The 
table also indicates that the assumptions about the magnitude of specific parameters 
underlying the present derivation are fulfilled, with the exception of Hq/k^ and Pe^ k. 
The latter is significantly larger than unity for the upper limiting case of Apo = lOOPa 
and K = 100 (/iq ~ lO/im). For this case, the advective fiux will distort the Boltzmann- 
distribution of the ions, although this is commonly neglected in most studies even for 
larger ho and Apo- For values of Apo < O(IO^) and 7c ^ O(IO^), Pen,k is indeed ^ 0(1). 
Furthermore, Pen fc > 1 refers to a regime where the Soret-effect and a temperature- 
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dependent permittivity have httle effect, as discussed before. Therefore, corresponding 
Hmitations are less important for the main conclusions drawn in this study. Given the 
definition of Ha, it becomes unbounded if Apo — > 0. In this case, instead of being 
the externally applied pressure difference, Apo needs to be replaced by the electrostatic 
pressure in the EDL, so that Ha = 1. The rescaling would not lead to the inclusion of the 



(now neglected) electro-migration terms on the LHS of (2.121 and ( 2.13 1, since they would 
still be 0{A^) smaller than the corresponding terms on the RHS. From this discussion it 
follows that values of Ha larger than unity are acceptable. Values of above 25 • 10~'^V 



formally contradict the DH-approximation. It is known ( Hunter||l993 ) and also verified 
in this work though that the DH-solution is surprisingly accurate even for much larger 
values of It fails primarily for k ^ 1 for which the EDLs of the walls opposite to each 
other overlap. The small values of S/E indicate that the modification of the streaming 
potential due to the Soret-effect is typically at least two order of magnitudes weaker 
than due to the permittivity change. For this reason, in what follows, these two effects 
are treated separately where necessary. 



For sufficiently large k (and thus also large k^*), tanh(K5-.) — )■ 1 and cosh ^(ks«) 



so that Fa 



k7' 



and Fp — )• 0. For k ^ 5, keeping the leading order term in Kg. 



> 
of 



each contribution, expression ( 2.34[ ) can be well approximated with 



st,p 



dxPo 



DH\K^5 
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1 ^liqCfs S*AT 



12 KS'DnT] Apo 



1 
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Ks* 



+ 



1 



2 KS'DnV 



(3.2) 



For the limiting case of infinitely small EDL-thicknesses compared to the channel height 
ho, it follows that k — )• oo and thus also ks* — >■ c». In this case, the ratio between the 
streaming potential at non-isothermal to isothermal conditions becomes simply 



isoth 



= fs.^exp{S*AT&) 



(3.3) 



DH\k- 



Integration in axial direction with a linear (dimensionless) temperature profile, O = X, 
subsequently leads to 



oth 



exp{S*AT) 
S*AT 



1 



1 + ^. (3^4) 

Hence, in this limit, the streaming potential increases approximately linear with in- 
creasing temperature-difference due to the Soret-effect while a temperature-dependent 
permittivity has no effect. In this case, the increase in streaming potential is due to the 
average increase of the EDL thickness with temperature by the factor f^}. Therefore, for 
very large k, it is sufficient to account for the Soret-effect by using an average increased 
EDL-thickness of (r-^)* = V[exp(S'* AT) - 1]/{S*AT). 

On the other hand, in the limit k — 0, one has Fa2 and F/^ — > 
one finds 



1/2. With (2.34) 



xH^st.,p 

dxPo 



(3.5) 

DH,-K^O 

This behavior agrees with the one obtained under isothermal conditions and is the result 
of two counteracting effects: decreasing k leads to a uniform charge density across the 
channel, potentially increasing the streaming current. Yet, k — > implies either /iq — 
or K — 0, where the latter commonly corresponds to Uoo ~> 0. The mass-averaged 
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velocity behaves approximately according to uq — ApQAho/{3riQ), which goes to zero for 
vanishing . Alternatively, for vanishing rioo , the charge density goes to zero so that the 



net result of k — is a vanishing streaming current (as can be verified with(2.22 1). Thus, 



at steady-state, the diffusion current expressed by (2.301 has to be zero as well and one 
can directly deduce (A0sf^p/Apo)K^o — > 0. 

Figure [2] depicts the ratio {A(j)st,p) / {A(j}st,p)isoth of non-isothermal/isothermal stream- 
ing potentials (from here on also called the 'streaming ratio') as a function of k for two 
temperature differences ((a),(b): AT = lOK; (c),(d): AT = 25K), two values of the Soret- 
coefficient ((a),(c): S* = IQ-^K-^] (b),(d): S* = lO-^R-^) and with AT/Apo and C\s 
as additional parameters. In (a) and (b), in agreement with the DH-approximation, the 
selected wall ^-potentials do not exceed C|s = 25 • 10~'^V (corresponding to C|s ~ 1); the 
results shown in (c) and (d) are obtained for much larger ^-potentials, i.e. 25 • 10~^V ^ 
^ 125 • 10^'^ V. In all cases shown, the dielectric permittivity is assumed to be con- 
stant {eT = 0). For estimates based on the DH-approximation (plotted with plain lines). 



numerical integration in axial direction of (2.34) was used to compute A^^f^p. Unlike the 



latter, {A<l)gt p)igoth is independent of the lateral coordinate and can be directly deter- 
mined by the same equation with AT = so that Kg. = k. Similarly, at higher and 
selected values of k, the streaming ratio was determined at by numerical integration in 
axial direction of ( |2.33 1 and marked with symbols in Figure [2] The pressure difference 



Apo is set to a constant absolute value of lOPa but changes sign in (a) and (b) . For pos- 
itive ratios AT/Apo, the thermal and mechanical gradients point in the same direction, 
while negative values denote the opposite. 

For the different values of the dimensionless parameter S*AT, the evolution of the 
streaming ratio with increasing k is qualitatively the same and can be divided into three 
distinct regions: the ratio approaches unity for small k, for which the EDL-thickness 
is very large compared to the channel height. Under this condition, a thermally in- 
duced modification of the already large EDL-thickness is of no importance and the non- 
isothermal streaming potential is identical to the isothermal one. For very large k, the 
pressure-induced streaming dominates and the ratio attains the value of [exp(S'*AT) — 
l]/(5'* AT). In this limit, as discussed before, the increase in streaming ratio is equivalent 
to the one obtained if the local, thermally induced expansion of the EDL is averaged along 
the channel length. In addition, the gain of A^^t p relative to {A(j>st,p)isoth is independent 
of the sign of AT/ Ap^ as well as of Ci^. As shown in the inset of Figure [21 (d), for large 
K and the parametric domain considered herein, Aip^t.p can be up to 16% higher than 



{A(j)st,p)isoth- For intermediate values 7t = 0(1), the Soret-term in (2.34) can be of similar 
order as the pressure- induced term and the value of the streaming ratio becomes sensitive 
to <^\s as well as to the sign and value of AT/ApQ. This general behavior is observable for 
all ^-potentials. At small values of C|s [Figure [2] (a) , (b)], the estimate based on the DH- 
approximation agrees very well with the full numerical solution except for k < 1 (shown 
is only C|s = 25 • 10~'^V but the same applies for smaller C|g). But even in this regime, 
the discrepancy is not larger than about 2%. At large values of C|s [Figure [2] (c), (d)] and 
large k > 5, the DH-solution is a still reasonable approximation of the (more accurate) 
full numerical solution. On the contrary, for small k < 1, the DH-solution erroneously 
locates the maximum in the streaming ratio at too large values of 7t. As discussed in sec- 
tion 2 and appendix [B] this is due to the overprediction of the EDL-potential-overlap by 
the DH-approximation. In addition, the magnitude of the maximum in streaming ratio is 
severely underpredicted by a factor of 2-3. At large ^-potentials, the DH-approximation 
generally tends to overpredict the streaming potential under isothermal conditions. This 
artificially lowers the relative increase in streaming ratio caused by the Soret effect. 
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As shown in Figure [2] (a) and (b), for small k, positive values of AT/Apo diminish the 
streaming potential. In this case, the pressure-induced advective ion transport is opposite 
to the ion transport by electro-migration due to the axial EDL-field generated by the 
thermally induced non-uniform expansion of the EDL-layer. It is important to stress that 
-counterintuitively- this opposing flow of ions is the reason for the reduction in streaming 
ratio rather than different thermo-diffusive fluxes of the ion species. As mentioned, both 
ion species are assumed herein to have the same coefficients of thermodiffusion, implying 
that the anion and cation number densities vary to the same extent in the direction of 
the temperature gradient. Thus, this effect does not alter the charge density and does 
not invoke any potential difference along the channel. Instead, the extension of the EDL 
vertical to the main flow direction varies in axial direction along the temperature gradient. 
Together with the strong variation of the EDL-potential within the layer, this causes a 
large axial field (within the EDL), which drives the ions in the EDL and in turn also the 
surrounding liquid. For negative AT/Apo, both ion transport mechanisms are additive 
and the streaming ratio increases. The maximum change of the streaming ratio increases 
from a little more than 1% at S*AT = 10^'^ and small C-potentials [Figure [2] (a)] up to 
more than 20 times at 5* AT = 2.5 • 10"^ and a large (-potential of C\s = 125 • IQ-^V 
[Figure I (d)]. 

3.2. Thermoelectric potential 



From (2.33), maintaining the thermal gradient while setting the pressure difference to 

(3.6) 



zero, one obtains the thermoelectric potential of a confined symmetric electrolyte: 

dx(f>st,T ksToo Ist.e dxQ 



Within the DH-approximation (for low (-potentials) , one finds with (2.34 1: 

%4'st,T 



1 ^iiqQ\s J J, tanh(Ks. ) tanh (K5. ) 




1 ^liqC\s 

2 D„Tj 



3ks' 

F„2 1 (3.7) 



The expression holds under the assumption that both ion species have the same ther- 
modiffusive mobility, i.e. that a diffusive charge separation process equivalent to the 
Seebeck-effect found in solids is absent, see appendix [Xj In the present case, the ther- 
mopotential is solely a confinement effect since it vanishes for 7t — > 00 (i.e. also Kg. — ^ 00). 
It also vanishes in the opposite limit, for k — > 0. 

Figure [3] illustrates A(j)gt,T as a function of k = k/iq; (a) and (b) show Acfi^t^T relative 
to the dimensionless Soret-parameter, 5'* AT, with ex = (i.e. E* — 0), while (c) and (d) 
depict A(f)st,T relative to the non-dimensional parameter E* AT, measuring the thermally 
induced change of dielectric permittivity. The Soret coefficient is set to zero in the latter 
case, i.e. S* — 0. In Figure |3] (a) and (c), the absolute value of the (-potential is varied 
from 5 to 25 • 10"^ V while it is varied from 25 to 125 • 10"^ V in (b) and (d). Solutions 
according to the DH-solution ( |3.7[ ) (numerically integrated in axial direction) are plotted 
with plain lines. Numerical solutions of ( |3.6[ ) (numerically integrated in axial direction) 
are displayed with symbols indicating each evaluation point. The Soret-coefficient used 
is S* — 10^'^K^^ ((a) and (b)) and the relative change of dielectric permittivity with 
temperature is E* = —5.1 • 10~'^K~"'^ ((c) and (d)), respectively. The temperature dif- 
ference is AT — 25K. Nonetheless, it was verified that the ratios A(pst,T / {S* AT) and 




Figure 2. Ratio {A(j>st,p)/{^4'st.p)isoth of non- isothermal/isothermal streaming potentials as 
a function of k = nho with AT/Apo and as parameters and with a constant permittivity, 
et = 0. In case of full numeri cal ev aluation ('Num.', lines with symbols), A</)st,p was computed 
by numerical integration of (2.331 in axial direction while the same equation was used with 
AT = (ks» = k) to directly calculate the constant value of {A<l)st,p)isoth ■ Each symbol 
marks the selected Tt-values for which the full numerical evaluation was performed. Estimates 
based on the Debye-Hiickel-approximation ('DH', plane lines without sy mbols ) were computed 
by numerical integration in axial direction of the analytical expression (2.341, divided by the 
corresponding isothermal value, (a) AT = lOK, S* = lO'^K"^ (b) AT = TUK, S* = 10~^K-\ 
(c) AT = 25K, S* = W'-'K"^ (d) AT = 25K, S* = 10~''K~\ inset: {A(l>st,p)/iA<l>st,phsoth 
obtained for 10 < k < 100 



A(l)st,T/{E*AT), are little affected by the value of 5* AT, respectively £;*AT, used for 
evaluation, at least in the parametric range considered herein. 

As can be seen in Figure [s] (a) and (b), as expected the thermoelectric potential based 
on the Soret-effect in confined geometries vanishes for very small and very large values 
of K. It reaches a maximum if half of the channel height, Hq, is of the same order as the 
EDL thickness e. at K « 1. This can be understood as foUovifs: the varying 

thickness of the EDL with temperature leads to an EDL-potential depending also on 
the axial coordinate x. While the EDL-field in axial direction remains constrained to the 
EDL itself, it propels fluid by electrohydrostatic and electroosmotic action. The afflhated 
convective charge transport subsequently leads to a measurable thermoelectric potential. 
Increasing ho beyond does not enhance the thermo-electro-osmotic fluid propulsion 
in the EDL due to the Soret effect. By contrast, the conduction current increases with 
increasing cross section, leading to a decrease of the thermoelectric potential. For the 
opposite limit ho <C the excess ion distribution across the slit channel becomes 

uniform and the local modification of the EDL-thickness by the Soret-effect is of no im- 
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portance so that the thermo-electro-osmotic fluid propulsion vanishes. As displayed in 
Figure [3] (a), for small ^-potentials, solutions according to the DH-approximation are 
practically indistinguishable from the full numerical solution. Notable is the relatively 
sharp peak in thermoelectric potential and its limitation to a relatively narrow K-range: 
for = 25 • 10~'^V, the thermoelectric potential l\4'st,T / {S* IST) at k = 10 is already 
only a little more than 10% of the maximum value at k = 2. In comparison, for much 
higher ^-potentials such as displayed in Figure [S] (b), the thermoelectric potential due to 
the Soret-effect is present within a much broader range of k. For these high (^-potentials, 
while predictions based on the DH-approximation deviate significantly from the full nu- 
merical solution for Tt < 4, still reasonable agreement is found for ^ 100 • 10~^V and 
K > 4. Both methods predict a strong change in curvature at k sa 4. With lirnax ~ 2 in- 
stead oiKrnax ~ 1-5 as obtained from the full numerical solutions, the DH-approximations 
slightly misplaces the location of the peak value. In addition, the DH-approximation over- 
predicts the maximum in thermoelectric potential by up to 40% (|C|sl = 125 • 10~'^V), 
and the maximum as predicted by the full numerical simulations is less pronounced. 

To quantify the effect of a temperature-dependent dielectric permittivity on the gen- 
eration of thermoelectric potential A(j)st,T: the latter is shown in Figure [3] (c) and (d), 
relative to the dimensionless parameter E*AT = eT^T/euq. The Soret-effect is absent 
{S* — 0). As discussed above, A(j>st,T/ {E* AT) vanishes for k ^ and for k > 5, inde- 
pendent of the ^-potential, leading to a relative maximum at k 1. As shown in Figure 
|3](c), for small ^-potentials the agreement between the full numerical solution and the 
DH-approximation is very good. For larger ^-potentials [Figure |3] (d)], the DH-based 
solution deviates significantly from the full numerical solution so that the peak value of 
A(j)st,T / {E* AT) is predicted to be up to (almost) twice as high as obtained from the 
numerical simulation (|C|s| — 125 • lO^'^V). In addition, the range of k corresponding to 
non-zero values of Acfigt^T / {E* AT) is more restricted around k « 1 than predicted by the 
numerical solution. This behavior is particularly visible at larger ^i^-values and similar 
to the one observed for A(j)st,T/{S* AT). Unlike the evolution of A<^st,T/ {S* AT) with k, 
remarkable is the fast approach of A(f>st,T/ {E* AT) to zero with increasing k > 5. While 
A4)st,T / [S* AT) maintains a non-zero value even for k = 100 (at least for larger 
A4>st,T I {E* AT)\-fi^^ ~ for any value of C,\g and independent from whether the solution 
was obtained from the quasi-analytical DH-approximation or from the full numerical 
model. 

To better understand the dependence of A(j)st,T/{E* AT) on k, it is instructive to 
consider the two contributions of a temperature-dependent permittivity (enq — f(T)) 
individually: Firstly, compared to isothermal conditions and as expressed in section [2.21 
Qiq = f{T) adds an additional term in the Korteweg-Helmholtz force (KHF), namely 



(V(/))^eTVT/2. Secondly, according to (2.17), an axial temperature-gradient acting 



on an electrolyte with euq = f{T) (i.e. E* ^ 0) varies the local EDL-thickness. While 
the first effect modifies only the streaming current (compared to isothermal electrokinctic 



flow), the latter affects both the streaming and the conduction current. Expression (2.17) 
suggests that a temperature-dependent permittivity modifies the local EDL-thickness 
qualitatively in the same fashion as the intrinsic Soret-effect. Hence, the contribution to 
the overall thermoelectric potential due to this effect alone should be qualitatively the 
same as observed for the Soret-effect. Indeed, artificially removing the e^-term from the 
KHF leads to similar results (not shown) as obtained for the Soret-effect (as illustrated in 
Figure [3] (a) and (b)). Consequently, the additional term in the KHF must be responsible 
for the sharp approach of A(j)st,T/ {E* AT) to zero for k > 5. This can be reasoned as 
follows. Within the DH-approximation one can show that the streaming current induced 
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in a non- isothermal channel by e = alone can be expressed by 
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(3.9) 



Here, I st.E* ,edl denotes the streaming current per non-dimensional thermal gradient 
9x0 due to the expansion of the EDL-thickness caused by e = f{T). Similarly, I st,E* ,ktf 
describes the streaming current due to the addit ional term in the KHF. The expression 
for Ist,E'.EDL corresponds to the third term of (2.25) if Js* is set t o unity and if S* is 



replaced by E* . Furthermore, Ist,E',KTF can be deduced from (2.25) if in addition dxPg 
and dx^st are set to zero and the /st, B'.BDL-term subtracted from the last term in ( 2.25 ) 



(which is the sum of Ist,E*.EDL and Ist,E',KTF)- For large k, Ist,E'.EDL approaches a 
constant value, while the sum of Ist,E',EDL and I st,E' ,khf approaches zero. Thus, for 
large k, the KHF leads to an opposing streaming current, which diminishes the streaming 
current induced by the expansion of the EDL due to e(T). 



4. Conclusions 

A non-isothermal, semi-analytical electrokinetic transport model of a symmetric dilute 
electrolytes in a slit-microchannel subjected to axial gradients in pressure and tempera- 
ture was developed. Unlike conventional studies of the Soret effect in bulk electrolytes, 
electroneutrality was not enforced, and it was assumed that the thermodiffusive mo- 
bilities of each ion species are the same. In the framework of the Debye-Hiickel (DH) 
approximation, an analytical expression for the flow field under the combined effects of 
axial pressure, temperature and electric potential gradients was derived. The derivation 
relies on a perturbation expansion in a small parameter, the height-to-length ratio of 
the microchannel. Solutions within the DH-approximation were compared to full nu- 
merical solutions. The temperature gradient was applied either parallel or anti-parallel 
to the main flow direction and parametric variations of the temperature- and pressure 
difference, ^-potential and of the intrinsic Soret-coefficient were undertaken. For diffusion- 
dominated ion transport and constant ^-potential along the channel walls, it was found 
that the electric double layer thickness (EDL) increases exponentially with tempera- 
ture. In comparison with isothermal conditions, this has two effects on the electrokinetic 
streaming potential: for large channel heights compared to the EDL-thickness, the on- 
average expansion of the EDL leads to an increase in streaming potential of up to 16%, 
independent of the orientation of the thermal gradient dxT with respect to the exter- 
nally applied pressure gradient dxPo- For channel heights of the same order as the EDL 
thickness, the ratio ('streaming ratio') between the non-isothermal streaming potential 
to the one obtained at constant temperature becomes dependent on the orientation of 
dxT relative to dxPo and, unlike for wider channels, sensitive to the specific value of the 
^-potential at the wall. For dxT/dxPo > (parallel), the streaming ratio is decreased 




Figure 3. Thermoelectric potential A(/)st_T of a confined sym metric electrolyte as a function of 
K = k/io and (^-potential (^g. DH-solutions accordi ng t o (3.71, integrated in axial direction, are 
displayed with plain lines; numerical evaluations of | |3.6[ ), integrated in axial direction, are plotted 
with symbols, where each symbol marks a point of evaluation, (a) A(l)st,T / {S* AT) (et ~ 0) for 
-C|. = [5,15,25] • IQ-^V, (b) A(j)st,T/iS*AT) (et = 0) for -C|s = [25,50,75,100,125] ■ IQ-^V, 
(c) A0,t,T/(S*AT) {S* = 0) for -C|. = [5,15,25] • lO'^V, (d) A<l>,t,T/{E* AT) {S' = 0) for 
-Qs = [25, 50, 75, 100, 125] • 10"^V. 



while it is increased up to more than twenty- fold for negative values (anti-parallel). This 
counter-intuitive result is a consequence of the non-orthogonality of the electric isopo- 
tential lines and the streamlines within the EDL. In other words, the thermally induced 
expansion of the EDL along the channel leads to an electric potential gradient within 
the EDL in axial direction, causing a thermo-electro-osmotic fluid propulsion. Depending 
on the sign of dxT/dxPo, this additional fluid motion either contributes or opposes the 
charge separation induced by pressure-induced electrokinetic streaming. For vanishing 
external pressure difference and small channel heights, the same thermo-electro-osmotic 
fluid motion causes a finite thermoelectric potential, although the thermodiffusive mobil- 
ities of each ion species are assumed to be identical. This distinguishes the present study 
from investigations of thermoelectric potentials in bulk electrolytes, which has a purely 
physico-chemical origin. The thermoelectric potential observed at confined conditions 
arises from the electro-hydrostatically and electro-osmotically induced ion streaming due 
to an electric field of thermal origin in the EDL. Finally, the consequences a temperature- 
dependent dielectric permittivity has on the streaming and thermoelectric potential was 
compared with those imposed by the intrinsic Soret-efFect. While the former is theoret- 
ically up to two orders of magnitude larger than the latter, in practice it only becomes 
important for very small channel heights (smaller than the Debye-length) , i.e. in cases 
with pronounced overlap of the electric double layers. At larger channel heights, the DH- 
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approximation achieves remarkable agreement with the full numerical solutions, even at 
large ^-potentials, but it suffers from inaccuracies in regimes with a strong EDL-overlap. 
The findings are useful for addressing the potential of small-scale fluidic waste-exergy 
recovery units, where the electrokinetic streaming is driven thermally. 



Appendix A. ThermodifFusion-Potential 

In the following, a brief overview of the conventional theoretical treatment of thermal 
diffusion potentials observed in multi-component bulk electrolytes subject to a tempera- 
ture difference is given. This summary is included in this study for the following reason: 
Starting point is classical nonequilibrium thermodynamics based on the phenomenolog- 



ical Onsager theory (de Groot & Mazur 1984), which is applicable to a wide range of 



different problems. As remarked by Hartung (2007), classical textbooks treat the Soret 



effect rather generally as one among many others, without going into great detail. Herein, 



a comprehensive outline of the matter is provided. As discussed in the book of Fitts 



(1962), it is assumed that the open system contains one ("second-law") heat flux q and 
P = K +1 material fluxes j^,. These vectorial fluxes are driven by P + 1 conjugate forces 
Xk, where Xq = Vln(T) and = V'^^Vi = 1, -P)- The latter denotes the spatial 
gradient of the chemical potential at isothermal conditions. The material fluxes are not 
independent of each other but obey jp — —J^k^i^k- electrochemical systems, it is 
useful to relate the material fluxes of the solutes to the motion of the solvent. For a 
system in mechanical equilibrium (i.e. the Gibbs-Duhem equation is valid), q and the K 



solute fluxes can be described by [page 66 in Fitts ( 1962 )] 
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L;,oVln(r)+^L^,;V(^)/i:, 



(Al) 



K 



-jfe = 4oVln(T) + ^4,V(^V:, 



(A 2) 



where L'^j are the phenomcnological coefficients. It is herein 



(A3) 



with V(/ii)|y — ^ fJ-i — (9T/ii)|pVT and fii is the chemical potential per unit mass 
of component i. The Faraday constant as the specific ion charge is denoted by F, Mi 
is the molar mass and (j) is the electric potential as before. With p^. as the mass of 
species k, ruk per total volume V, one has fii = ^ii{T,p, pi, ...px) so that ^{pi)\T = 

dp{pi)\T,pi'^P + Y^iLi dpiipi)\T,p^Pi with dp{pi)\T,pi = ^ ttie partial volume Vi per 
mass rrii. Hence, 



K 



j,=4oVln(r) + ^L',, 



V,Vp + Y,dpM)\T,p'^Pi + 

1=1 * 



(A 4) 



.p, while L'fco 



The diffusion coefficients are defined with Dn,ki = J2i=i ^ki'^piiPi)\T,i 
pkTDT,k with DT,k as the thermodiffusion coefficient of component k [page 79 and 102 
in Fitts (1962). Note that, compared to the latter, in this summary the order of the 



indices is k,i,l instead of i,l,k]. Neglecting pressure-induced diffusion and defining ion 
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mobilities ujk with evkPk^k = X^iLi ^'ki^i^ I ^'U leads to 



K 



jfe = pkDT,k^T + ^ Dn,ki'^Pi + evkPkUJk'^4'- 
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(A 5) 



1=1 



In the present case K = 2. From (A5), expression (2.1) can be found by neglecting cross 
diffusion due to concentration gradients, approximating cj) k, ip (Osterle) and inserting 
the result into dtUk = —NA/MkV ■ j^, where Na is the Avogadro-constant. Note that 
pk = nkMk/NA and F — cNa- The intrinsic Soret-coefhcients are defined by pkF>T,k = 
Ylf=i PiSiDnM so that 



K 



jfc = ^ Dn,ki (Vp/ + piS'iVT) + evkPki^k'^ (i 



(A( 



1=1 



Alternatively, each material flux expressed with ( A 2 1 is weighted with parameters 
to be determined, summed over all K and the result subtracted from (Al). This leads 
to 

9 = E ^^Jfe - ( ^00 - E ^l^'f^A ^li^(^) - E ( ^« - E ^^^'h ) "^'"^V:- (A 7) 

k=l \ k=l / 2=1 \ fe=l / 

The weighting factors Ql, the so-called heat of transports, are now selected so that the 
last term in (A 7) vanishes, i.e. 

K 



^Oi — E ^k^'ki 



fc=l 



Defining the overall heat conductivity A via 



K 



XT = L'r 



00 



feO 



fe=l 



leads to 



K 



q = Y.Ql3k-^^T. 



(A 8) 



(A 9) 



(A 10) 



fc=i 



With (A8| and the Onsager-relation L^^, = LJ,- so that Lq^ — L'^q, one has L'^ 



S^i-^fciQ*- Consequently, with (A3), the material fluxes can be expressed by 

K 



kO 



■3k 



i=l 



Q*VHT) + V {p. 



(All) 



In solution chemistry, it is a common practice to express the dependence of the chem- 
ical potential on the composition in terms of the activities = "fifhi, with 7^ as the 
activity coefficients and rhi — Ui/psoiv as the molarities. One has V(/ii)|T = V^Vp -f 
i?T/MiVln(ai), where R is the ideal gas constant. This expression omits the depen- 
dence of the chemical potential of species i from the concentration of species k, i.e. 
cross-diffusional effects are implicitly removed. This is reasonably accurate for very low 
concentrations of the solutes. One obtains ( Hills et aT]|1957 ) 



K 



' Jk — E ^'^i 



(A 12) 
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Neglecting the pressure-induced diffusion and assuming constant activity coefficients (di- 
lute limit), one finds with R = NAks and Q* = NaQ*/M.i 



3k 



K 



1^1 



ki 



ksTN, 



A 



VT + Vn, 



'ksT 



(A 13) 



As mentioned, this equation can only be derived by neglecting cross-diffusion between so- 
lute species, i.e. L^.^ = for k ^ i. This, according to the definition of Dn ^i introduced be- 
fore, yields D„,fc = D^^kk = L'f.^.dp^{fik)\T,p with {fik)\T,p = ^AksT /MtM{-iknk/ Psoiv)- 
Using dp^{.) = N A IMA J.) one finds dnAuik)\T-n = NlknT I {ukMi) and 



NA/Mkdn^{-) one finds dp^{pLk)\T,p 
Mk " 



-jk= D„ 



''Na 



ksT^ 



N\kBT/{nkMi 



ksT 



(A 14) 



In the absence of cross diffusion between different ion species and equivalent to the 
treatment chosen in the main part of this work, (A 14 1 provides the non-isothermal 
diffusive ion flux in terms of the (ionic) heats of transport (Wiirger 2010). From this 
equation one can also deduce that — Q^/ {kgT'^). Within the general framework of 
linear response theory for non-equilibrium processes, the particle mobilities should be 
independent of the thermodynamic conjugate forces, including T. In this case, for a 
non-isothermal system, (A 14 1 should read 

Mk \ Ql VT 



3k 



■Na 



nk 



knTnc, T 



Vrik + Uk 



Vke 

kuTnr 



(A 15) 



Linearizing the thermal conjugate force VT/T leads to S^. ~ Ql/ (ksT^). 

For a (mass-) closed system at steady-state {t — > oo), each material flux vanishes, i.e. 
jk = for all k, and thus VkFj k/ {Dn,kMk) = as well. Summing over all K leads with 
( |XT4| to 

K 



K 

E 

fc=i 



evkTik 



ksT 



VT 



evk_ 
kuT 



(V<i 



)\t- 



V^(ez/fenfe) =0 



(A 16) 



k=l 



{Vc, 



j|t— >-oo 



(A 17) 



The last term vanishes at electroneutral conditions present in the bulk electrolyte outside 
the EDL. The corresponding steady- state thermoelectric potential reads 

Y,k=i e-VknkQl (YZ\ 

From this equation one can deduce that the bulk thermoelectric diffusion potential, 
equivalent to the Seebeck-effect observed in metals and semi-conductors, should vanish 
for symmetric, dilute electrolytes, if the (ionic) heats of transport (or equivalently, the 
intrinsic Soret-coefficients) of each ion species are identical. 

=1 Jfe = 0. Defining the average salinity 



For completeness, jk — for all k implies 



with 2no = X^s^i "i: 

K 



/VT 



one finds with (|A 14|) 

Ql 



Eti^nk 



• k=l 



no 2kBT 



no 



K 

nk 



evk 



' j-^ no 2kBT 



= (A 18) 



Charge neutrality leads then to the salinity gradient at steady-state (Soret equilibrium) 
( |Wurger|[2010l ): 



no 



(A 19) 
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Appendix B. EDL-potential 



In this section, for completeness, the fuU analytical solution of expression (2.7) is 



repeated, using the nomenclature as pertinent to this study. Integrating (2.7) once in 



vertical direction (K5. is a function of X only), using the symmetry condition and ^E* 
at Z = 0, leads with cosh(*) = 2cosh^(^/2) - 1 and Z = 7tg»Z to 



±[cosh^(^'/2) -cosh-'(«'|c/2)]i/2. Following Burg reen & Nakache (^1964L one substitutes 
sin(7?) = cosh('I'|(,/2)/cosh(\I'/2), b — [cosh(^|c/2)J~^ and integrates again. This provides 

" (Bl) 



Ks'Z — b 

Jl-b^sm^i-d) 

where = 7r/2 = at Z = 0. Hence ( [Levine et aL]|1975[ |Keh fc Tseilil|2001[ ), 



Fi^,b) 



(B2) 



where F{'d,b) = d-d/yl — 5-^sin^(??) is the incomplete elliptical integral of the first 
kind ( Langmuir|1938 ). The local EDL-potential is given by ^ = 2arcosh[ cosh(\E'|c/2)/sin('i9) 
and is determined by the value of the elliptical integral at the wall {Z = 1), namely 
= F{Tr/2,b) ~Ks-/b so that ^-i^ = 2arcosh[cosh(C| J2)sin(i9|^)]. 
Within the present context, it is irrelevant from which variables the local Debye- 
parameter kq* is a function of. Hence, for the following discussion, ks' is substituted by 
K. In Figure 4I the distribution of the EDL-potential = '0(^)/C|s along the ver- 

tical Z-coordinate is depicted for two values of the ^-potentials, namely = 25 • 10~'^V 
in (a), corresponding to C|s ~ 1, and (\s = 125 • 10~'^V in (b), corresponding to ~ 5. 
For K = (0.5, 1, 5, 10), the exact solution (B 2 ) is shown in comparison with the numer- 



ical solution of expression (2.7) as well as with the Debye-Hiickel (DH)-approximation 
(2.9). Note that 4'£)//(Z)/C|g is independent of the wall ^-potential. The boundary value 

problem described by (2.7), along with ^'(±1) — (^^ at the walls, was solved by collo- 
cation with the BVP4(j-tunction implemented in Matlab (Version 7.11.0, R2010b) on 
a Dell Precision T7500 workstation operated with CentOS 5.8. A computational mesh 
with approximately 500 grid points in Z-direction was sufficient to obtain numerical 
solutions which match the exact solution up to an accuracy of 0.1%. The agreement 
between the exact and the DH-solution is good for = 25 • 10~'^V and larger values 
of 7c, as can be expected from the simplification made in the DH-approximation. For 
overlapping EDLs, i.e. for 7c ^ 1, noticeable disagreements occur, as particularly visible 
for = 125 • 10~^V. Furthermore, it is noted that the DH-approximation generally 
overpredicts the EDL overlap. 



Appendix C. Verification of self-propulsion of EDL 

In the main part of this work it is shown that, for a confined system, a finite ther- 
moelectric potential emerges, even if the thermodiffusive mobilities of the ion species 
are identical. Essentially, it is caused by an axial gradient of the EDL potential which 
propels the fluid by means of the electrohydrostatic pressure and the electro-migration 
force (Maxwell stress). For an isothermal elect rokinetic system, it is well-known that the 
EDL itself does not set the fluid into motion. It was shown by Levich (1962), page 484, 
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Figure 4. EDL potential = for k = (0.5, 1, 5, 10). Exact solution (EX,black 

solid line s) (|B 2[ ) in comparison with numerical solution (NM, thick dashed red lines with sym- 
bols) of |2.7| and Debye-Hiickel (DH, thin dashed green lines with symbols) approximation 
(a) C|« = 25 • 10-^V (C|, « 1), (b) C|s = 125 • lO'^V (C,, « 5). 



that, for such systems, the electric body force due to the EDL and the osmotic pressure 
contribution exactly cancel each other and that only the externally applied electric field 
is relevant for the fluid propulsion ( Pascall fc Squires]|2011 ). Therefore, the phenomenon 



described in this work may be suspected to be an artifact of the lubrication approxima- 
tion (LA) used in the analysis. In the following, without relying on the LA, it is shown 
that a non-isothermal EDL can indeed propel the fluid (while an isothermal EDL can 
not). 

For an incompressible Newtonian fluid with constant dielectric permittivity, the Stokes 
equation reads 

= -Vp + ?7oV2v-p/V(/). (CI) 



To eliminate the pressure, the curl-operator V x (.) is applied to ( C 1 ) , so that V x Vp = 0, 
and n = V X t> is the vorticity. If only the EDL potential ?A is present, then (j) = ip. One 
finds 

T^oV^n^ V X (p/W). (C2) 

The charge density p/ is a scalar, so that the product rule for the curl-operator and 
V X VV- = yields 

TyoV^n = (Vp/) X W. (C3) 



Replacing the charge density with expression (2.4 1 and evaluating the derivatives leads 
to 

7?o V^n = -2kBT^n^ [fg, cosh(«') V* x V* + sinh(*) ( V/|. ) x V^-] . (C 4) 
Since x V* = 0, one has 



, ^ -V^n = -2sinh(vl/) (V/|.) X V^. (C 5) 

In the present study, fs* is a function of the axial coordinate so that 7^9/ (?^oofcBToo)V^n = 
— 2sinh(^')9j. (/|.) d^'^ 7^ 0. On the contrary, under isothermal conditions, one has 
fs* = 1, leading to 77o/(^oo^s7oo)V^n = 0. This implies that in this case the flow 
field is unaffected by ^ for any value of 4' and, consequently, the EDL potential does not 
propel the fluid. 

It is emphasized that, in the present derivations and unlike in other studies of elec- 
trokinetics, an expression for the osmotic pressure of the ion cloud is never explicitly 
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added to the governing equations as a body force, neither in the derivation based on 



the LA nor in the one leading to (C5|. In the LA, a corresponding term automati- 



cally emerges from the momentum equation in vertical direction, (2.13). Its combination 



with the electromotive force in the EDL gives, under non-isothermal conditions, a non- 
vanishing force contribution. An equivalent force is obtained in the vorticity formulation. 
As shown above, the appearance of the propulsion due to the EDL itself can be traced 
back to the circumstance that, herein, the EDL-Maxwell stresses are not irrotational, 
while they appear to be under isothermal conditions. It needs to be addressed whether 
p f always has to be of a form so that the curl of the EDL-Maxwell stresses vanishes also 
under non-isothermal conditions. In other words, the result of (C5) might be an artifact 



introduced by a thermodynamically inconsistent expression for the charge density. As 
outlined in the following, the present form of simply implies that the ion cloud is not 
in mechanical equilibrium, while it is at constant temperature. 



According to Fitts (1962), page 43, the incompressible Newtonian Navier-Stokes equa- 



tion can be written as 



dtv - V^t; = -Vp + pqY, 



(C6) 



where po is the fluid density and Y is the sum of all body forces. The RHS of (C 6 1 can 



be expressed by the general form of the Gibbs-Duhcni equation for a multicomponent 



system (Fitts (1962), page 44) 



Vp - pY 



p 

E 

fe=i 



Pk 



(C7) 



The mass of component k per total Volume V is denoted by p^. The sum contains all 
components P = K + where K is the number of solutes. The spatial gradient at 
constant temperature of the overall chemical potential is defined according to 



(C8) 



where Hk — p'j^ — VkF/Mkil^ {Mk is the molar mass, F = cNa is the Faraday-constant with 
A'^ as the Avogadro-number). It is assumed that the gradient of the chemical potential 
of the solvent vanishes so that the effective upper limit of the sum in (C7| is K. As 



deduced before from (2.2), the thermo-electro-chemical potential of the ion cloud should 
be suitably described by /i^ (see plj ), i.e. p'j, = UkF/Mkip + kBTaoNA/Mk[Mnk) + SlT]. 
Thus 



ksTooNA Vnfc 



(C9) 



Employing Uk — Na/Mj-pu and "Y^^^i evk^k = X^i^i ^^^k^k — Pf (for the neutral solvent, 
'^p — 0) yields 



K 

E 

k=l 



Pk 



(CIO) 



with n = X^A^i "-fc- The last term on the RHS is the gradient in osmotic pressure, i.e. 



fcsTooVn. To simplify matters, for the further analysis low values of the ^- 



potential will be assumed. Under isothermal conditions, pf 
K are constants) 



-eiiqK^i}}, yielding {euq and 



(j^PkV^^'A 

\k=l J , 



(Cll) 



?oth 
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where N = n/rioo — 2cosh(\E') « 2 + VP^. Hence, at constant temperature, the RHS of 
(C 11) is identical to zero ( Levich|[l962 Squires &: Bazant^^2004j ), and the ion cloud is in 
mechanical equilibrium. By contrast, under non- isothermal conditions the ion cloud is 
not necessarily in mechanical equilibrium, initiating a propulsion of the fluid. 
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